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Abstract 

We investigate the gravitational radiation emitted by an isolated system for gravity theories 

with Lagrange density f{R) = R + aB? . As a formal result we obtain leading order corrections to 
O ■ 

\ the quadrupole formula in General Relativity. We make use of the analogy of f{R) theories with 

^ ^ ' scalar-tensor theories, which in contrast to General Relativity feature an additional scalar degree 

. of freedom. Unlike General Relativity, where the leading order gravitational radiation is produced 

by quadrupole moments, the additional degree of freedom predicts gravitational radiation of all 
O^l multipoles, in particular monopoles and dipoles, as this is the case for the most alternative gravity 

theories known today. An application to a hypothetical binary pulsar moving in a circular orbit 



o 
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yields the rough limit a < 1.7 • 10^^ m^ by constraining the dipole power to account at most for 
of the quadrupole power as predicted by General Relativity. 
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I. INTRODUCTION 



One of the most impressive endorsements of General Relativity Theory (GR) is the agree- 
ment of the predictions of the famous quadrupole formula for gravitational radiation with 
indirect measurements of the energy loss of binary pulsars. It is thus natural to test mod- 
ified gravity theories by deriving the corrections to the quadrupole formula and comparing 
them with experimental data. For many types of theories this has been done in the past 
Though this problem is still open for metric f{R) theories with an action 

S = J f{R)V^ d'x + 5m, (1) 

where in contrast to GR the Einstein-Hilbert Lagrangian density is replaced by a nonlinear 
function f{R). We will always assume the function / to be smooth. Sm is the standard 
matter action. In the past years, this type of theories has become very popular to heuris- 
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y gain insight in the problem of dark energy. For an overview one may consult e. g. 
3| and references therein. It is also worth mentioning that quadratic corrections aR^ to 
the Einstein-Hilbert action have been considered in the context of inflationary cosmology 
already in the 1980's jl, 5|. 

In this article we prepare the way to investigate the energy emission of binary systems by 
gravitational radiation. The basic equations of f{R) gravity are given in Section II. For our 
purposes it will be convenient to work in the scalar tensor formulation of quadratic f{R) 
gravity. In Section III we employ the linearised field equations of quadratic f{R) gravity 
to derive the weak gravitational fields emitted by a localised source and expand them into 
multipoles. The linearised f{R) gravity has been investigated for example in [2I, gHq]. For 
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a treatment of linearised scalar tensor theories see [ll, llOj • In Section IV we dwell on the 
energy-momentum complex in quadratic f{R) gravity as an analogue to the Landau-Lifshitz 
complex in GR. The leading order correction to the quadrupole formula in terms of momenta 
of the energy-momentum tensor is derived in in Section V. In Section VI we finally illustrate 
the correction with an application to binary systems in circular orbits. 

Notational conventions: Greek letters denote space time indices and range from to 3, 
whereas Latin letters denote space indices and range from 1 to 3. We take the sum over 
repeated indices within a term. 
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II. THE FIELD EQUATIONS 



Consider a 4-dimensional pseudo Riemannian manifold with metric gfj,u of signature 
(—,+,+,+). We write g = detg^^i, and denote the Ricci tensor of g^^, by R^^. The variation 
of the action ([T]) with respect to the metric yields the Euler-Lagrange equations 

f{R)R,. - ^-f{R)g,. - V,VJ'{R) + g,.UJ\R) = ^T^,, (2) 

where R = g'^'^R^^, T^u = {—'^c/^/^){6SM/Sg^'^) is the energy-momentum tensor, c the 
vacuum speed of light, G Newton's constant, the covariant derivative for g^^^, and Dg = 
V^V^. Taking the trace of (|2]) we obtain 

SDJ'iR) + f{R)R - 2f{R) = ^T, (3) 



where T is the trace of T^^. We now assume 

f{R) = R + aR^ (4) 

and make use of the equivalence between f{R) gravity and scalar tensor theory by defining 
the scalar field (p := f'{R)- This identification is feasible since f"{R) 7^ holds for our choice 
of f{R), and f'{R) is thus invertible. We define the scalar field (f hj (p = 1 + 2a(f, where 
we have chosen the asymptotic value such that a renormalisation of Newton's constant is 
redundant. Then the equations ([2]) and ([3]) are equivalent to 



(5) 
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= jT + —ip. (6) 

The field ip thus has the effective mass h/{c^/6a). Since we aim at the investigation of 
isolated systems at the scale of the Solar System, we expect the theory derived from (jll) to 
be a small perturbation of General Relativity. Thus the dimensionless quantity aR should be 
small compared to 1. This still holds if the parameter a varies appropriately with the scalar 
curvature and therefore the local matter density. For our quadratic model this fact reflects 
the concept of the chameleon effect pjj] , which states the possibility that the Compton wave 
length A = a/Go of the field (p is smaller or larger in regions with higher or lower matter 
density, respectively. In our specific case however, we assume the Compton wave length to 
be constant. In this sense theory should be viewed as a local effective field theory which is 
valid for a certain range of the parameters. 



III. GRAVITATIONAL RADIATION IN f{R) GRAVITY 



Consider weak perturbations of the Minkowski spacetime metric rj^i, = diag(— 1, 1, 1, 1). 
The metric can be written as 

9fiu = Vfiu + hf,^, (7) 

where the coefficients of the perturbation satisfy \hf^i^\ ^ 1. In what follows the indices are 
raised and lowered by r]^,^. For the field we have already chosen the asymptotic value 1, of 
which 2a(p is the perturbation. Moreover, the field equation (|6]) is inhomogeneous linear in 
(fi, so that the linearisation in the perturbations is simply achieved through the replacement 
of Dg by Drj. Let h = h^, define 



and choose the gauge 



7, 



0. 



Up to linear order in h^y and y?, equation (j5]) can then be written as 



-T 



(8) 



(9) 



(10) 



In the slow motion approximation and at large distances from the localised sources, a special 
solution of (fTOj) can be derived in analogy with the GR case. Since T^j, is divergence free, we 
can to express the spatial components of 7^1, in terms of the quadrupole momenta of T*^". 
Thus we obtain 



For the field ip we write equation (Q as 



d^x'T^%t-\x\/c,x')x"x 



(11) 



S, 



(12) 



where a := l/\/6a and 



S = T 



1 + 



+ 
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3a' 



-{VVf + UV + VW 



(13) 



is the source T (to leading order) extended by the terms which are quadratic in the pertur- 
bations. These are expressed in terms of the Newtonian and post Newtonian potentials U, 
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W and V, which are given by 12|, Il3 | 



U{x,t) 



W{x,t) 



V{x,t) 



AG /•(-2)r"0(a;',t)^3^, 



Sc^ J \x — x\ 
Sc^ J \x — x'\ 

' ; ; ax . 



(14) 



(^l _ e-"l— 'I j rfV, 



\x — X 



Here 2)2^00 jg ^^le leading order time-time component of T^'^. 

A special solution of f[T^ is the convolution of the source with the Green's function of 

n 

the Klein-Gordon equation [14| . 



6{t- \x\/c) "^1 (acVt^ - (|a;|/c)2j 6 {t - \x\/c 



\x\ 



(15) 



V^^-d^l/c)^ 

where Ji the Bessel function of first order and the Dirac, and Heaviside, distribution kernel 
is denoted by 6, and 9, respectively. Then 

2Ga^ 



^p(t, x) 



,3 ,S {t — \x — x'\/c, x') 



(TX 



\x — X 



(16) 



t—\x—x'\/c 



dt' / d^x' 



aJi (ac^/{t-t'y - {\x-x'\/c)^ 



^{t-t'y-i\x-x'\/cy 



S{t,x') 



If the source emits a single pulse, the field ip observed at a distance consists of this pulse 
diminished by the factor l/\x\ given by the first term on the right hand side of f|T8|) . and a 
wake represented by the second term. The scalar mode is thus dispersive; for a discussion 



of the dispersion of plane scalar waves see 
After the substitution 



t'=t 



\X — X 



1 + 



a^\x — X 



f\2 ' 



dt' 



a^c\x — X 



ds 



1 + 



;i7) 



in the second term of the right hand side of (fT6|) . the solution can be written as 

S{t-\x-xyi + j£^/c,x' 



ip{t,x) 



G 
3ac^ 



ds[Ji{s)e{s) -6{s)] / d'x 



\x — £c'L /l + 



:i8) 



|a;-a;'|2 

We assume the source S* to be a smooth function of time. In the limit a — )■ 0, the spatial 
integral in (fTSl) is independent of s, and because of Ji{s)ds = 1 we have lima^-o ^V^ = 0. 



Let r := \x\ and n := x/r. At large distances r from isolated and slowly moving sources, 
we can expand f lTSj) into multipoles. Since the scalar contribution to the energy flux (equation 
below) is quadratic in ip and 7^1,, and moreover, the leading order radiation predicted 
by GR is produced by quadrupole moments, we need an expansion for taking into account 
up to hexadecapole moments. We therefore derive the fourth order Taylor polynom of the 
integrand in (fT8|) around the origin in the variable a;', while ignoring the explicit dependence 
of on as'. Consider a source with an extension characterised by a typical length d <^r and 
moving slowly at a velocity characterised by a typical frequency uj <^ c/d. The polynomial 
can be written as a sum of terms of approximate order O ^^J^^r^sj = O ^^^j , m + l = n. 
In fact the expansion coefficients of the monomials x'^^ ■ ■ ■ a;'*" are not proportional to r"™, 
but depend on r also by means of the function 



For convenience, we will ignore the function p when making use of the Landau Symbol O. 
We define the retarded time 



and the distribution kernel 



T{t,s):=t-—- (20) 

P{S)C 



q{s) := p{s)[Ms)eis) - 6{s)]. (21) 



We are then left with 



cp{t,x) = [ dsq{s) [ d^x'\l + Fi{s)x" + Fij{s)x"x'' (22) 

SiT{t,s),x') + 



+ Fijk{s)x"x'^x"' + 



d', 
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(n + m = 5), where 
FAs) 



Hi 



r 



c dt 



(23) 



Fijkis) 



riirij 



2r2 ^ rc 9t 2c2 



5ij 



2^3 2r2c dt 

^3 



+ 



p2(s) p(s) 9 

^3 



2rc at 



2rc2 6c3 



3p (s) ^ p (s) d ^ p (s) & 



Fijki{s) 



niUjUkni 

- UiTljdkl 

— ^ij^kl 



2^3 2r'^c dt r(? dt'^ 
35pS(s) 35p^(s) d lbp^{s) 



p^{s) 



d^ p\s) d^ 



8r^c dt 



+ 



15p6(s) ^ 15p5(s) 9 ^ 3p^{s) 92 



8r2c2 12rc3at3 2Ac^ dt^ 

p^{s) d^ 



4^4 



4r3c at 2r2c2 5^2 12rc3 dt^ 



3p\s) 3j^d_ p\s) d^ 



^r** Sr^c at 8r2c2 at2 
In what follows we drop the post Newtonian source terms that are quadratic in the pertur- 
bation fields. Note that these contributions would lead to corrections in the energy emission 
formula which are quadratic in G. Therefore we will neglect them and proceed with T^^ as 
the main contribution to the source, i.e. S = rj^^T'^''. If we consider a perfect, non viscous 
fiuid with mass density p, pressure V and velocity field v = {v^,v'^,v^), we have 



[p{t,x) + 0{c~')], 
c[p{t,x)v' it,x) + 0{c~^)] 



T°° (t, x) 
(t, a;) 

(t, x) = p (t, a;) v' (t, x) (t, x)+V (t, x) 6ij + O (c^^) 



(24) 



We express the spatial integrals over the source in (JTTi) and (1221) using the following momenta 
of the energy-momentum tensor, 
1 



M'"(t) := 



and the quantities 



/ d^xT^° {t,x)x^", S'^^-{t):= I d^xT'^t,x)x^" 



(25) 



(26) 



Here denotes a string of indices . . . n = 0, 1, 2 . . ., and x^" := 



X 



X*". Moreover, 



we will denote := A^ooo- It is useful to introduce the following linear combinations of 
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the quantities (l26i) . 

D\t) := M',M+M\,,{t), (27) 

QYit) ■■= lM%,{t) + MUt) + lMg,{t), 

Ofis) ■■= lM'ioit) + lMg\it) + lMf,it) + lMf,{t) 
Of{s) := - Sf,{t) - Sf,{t) - \Mt,{t) - - Mf,{t), 

QC, QC, IK 1 1 

Hf\s) := - - S^^it) - Sg^^t) 

Taking into account up to quadrupoles, and hexadecapoles, for 7*-', and respectively, the 
asymptotic fields ffTTl) and fl22|) can be written in terms of the quantities fl25l) . fl26|) and fl27|) 

7^^(t,^) = ^^^^^(^--Z^)' (28) 



G 



M (t) + m D' it) + rurij Q'l {t) + 5^^ Q'i (t) (29) 

where we have used the same symbols 7*-' and (p for the approximations. From fl2^ we 
infer that the dimensionless field ay? depends on a through the function p. This irrational 
dependence on r arises by the same reason as the Yukawa like terms e"*"/^^ in the solutions 
for (p for an isolated system, since the field (p has a range 1/ \/6a as per equation ([6]). In 
the same way as in the 1/c expansion of the metric field 13], the approximation is 
parametrised by a dimensionless parameter a/£^, where £ is a typical length depending on 
the model. As for the 1/c expansion of the field of a quasi static isolated source, also in 
the present case I corresponds to the distance from the source r, i.e. the distance the wave 
has propagated through. Consequently, even the field ray? converges to zero for r — cxd, 
that is, other than the pure metric field (128|) . the scalar field decays faster than 1/r. This 
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is expected since the scalar field is massive, leading to a nontrivial dispersion relation for a 
propagating wave. 



IV. THE ENERGY MOMENTUM COMPLEX 



In order to derive the energy flux of a gravitational wave in f{R) inspired scalar tensor 



theory, we need an analogue to the Landau-Lifshitz complex t'^j^ in GR 15|. This can be 
obtained by using the method in [l6^, where an energy- momentum complex in the Brans- 
Dicke theory is presented. This method applies to theories derived from convex or concave 
functions /, i.e. theories that are equivalent to scalar-tensor theories. In this section, we 
consider a general metric tensor g^^i, and scalar field (p. In particular, we do not constrain 
the asymptotic values of these fields. For an alternative derivation of the energy-momentum 
complex cf. j^. 



17|. 



Defining 



and 



we can write the generalisation of the Landau-Lifshitz complex as 



SttG 



-9) 



(30) 
(31) 

(32) 



By construction, the energy-momentum conservation laws then can be cast into the form 

d, {<P{-g) (T'^^ + tn] = 0. (33) 
Using fl30|) and fl3T]) , it is straight forward to express the energy momentum complex f l32|) in 



terms of the fields g^^i^, 0, their first partial derivatives and the connection coefficients F 

as 



A 

fJ,U 



28^4)8''^ 



1 



2 /^1/ 



(34) 



+ 



SvrG 



+ (9>FL + 9>F^J g^'^ - {d^<Pg''^ + d''<Pg^^) F^, 
- {g^^T\^ + g'^^Tl^) 9^ 
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The energy flux in an arbitrary direction is given by the component t' 



V. ENERGY EMISSION OF ISOLATED SYSTEMS 



Consider a plane gravitational wave propagating in vacuum in the direction. In addition 
to the gauge iQ, it is possible to perform a further gauge transformation which makes the 
(p independent part of the perturbation transverse and traceless (TT) 2, 6, oj, such that we 
can write 

h'^'it - x'/c) = 7^^^(t - x'/c) - 2ar]f"'<f{t, x') (35) 

In this gauge, we evaluate the energy flux to leading order in the perturbation flelds ■jf^u and 
if. By angle brackets we denote the average over a four- dimensional spacetime region with 
an extension which is much larger than a typical wavelength. The formula then yields 



Oi 



(36) 



327rG 

The flrst term on the right hand side of fl36|) can be evaluated in the same way as in GR by 
means of (128|) and the trace-free quadrupole tensor 



1 



(37) 



From f l29|) we obtain for the second term on the right hand side of fl36|l the asymptotic value 

8G^ 



24:a\doipy 



3c^r 



+ 2niMD' + n,nj (2MQ^i + D'D^j + 25ijMQ^ 
+ niu.nuni (2MHf' + 2D'df + Q'^Qf 



(38) 



+ 5,,5ui (2MHf' + Q^iQ 



The total power of the source is obtained by integrating ( 136|) over a sphere with radius r: 



3^^ +3 



M [ 2QX + GQ'i ) + D'D^ 



(39) 



+ 



15 L 



+ gf gf - logi' ) + 2gi-'gi^ + i5g^2^g^/ 
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where the dot denotes the derivative with respect to t. Note that the additional degree 
of freedom represented by the scalar field predicts radiation of all multipoles, in particular 
monopoles and dipoles. These are of lower order than the original quadrupole contribution 
in GR represented by the moments Q*-', and could principally lead to a non negligible 
contribution in concrete applications. Since the quantities M.^ and as well as their 
derivatives vanish in the limit r — i- oo, the energy radiated by means of the scalar degree 
of freedom is absorbed as the field (p disperses completely when propagating to timelike 
infinity. 



VI. APPLICATIONS 



The most interesting application of the formal result fl39|) is the energy loss of binary 
systems by the emission of gravitational radiation, in particular binary pulsars such as the 
PSR J0737-3039 system [JJ, 



\\\ . In quadratic f{R) gravity, the non-relativistic motion 
of compact objects is governed by the Newtonian potential with an additional Yukawa 
correction. This implies that the Keplerian orbits also need appropriate corrections, cf. for 
example 12|. A general treatment including generic orbits is beyond the scope of this work. 
In order to obtain a rough estimate for the correction terms in (!39|) . we apply the formula 
to the radiation of binaries in circular orbits. For this application we expect that the total 
mass of the system changes on a time scale which is much larger than the orbital period. 
Hence, the monopole contribution to ( 139|) is negligible, ~ 0, whereas for appropriate 
mass densities the dipole term is the leading order contribution. 

Consider a binary system consisting of two point masses mi and m2 in a circular orbit 
moving at angular velocity u. For mi ^ m2, the dipole moment does not vanish. The only 
non-relativistic correction is a modification of Kepler's third law 12| . 



CO 



G{mi + 7X12) 



1 + 



d 



'6a 



-d/ y/Oa 



(40) 



which can be numerically solved for the orbital separation d. This non-relativistic correction 
is parametrised by a/rf^, while the corrections arising from the propagation of the field are 
a function of a/r^. The choice of the ratio d/r is constrained to be small for the multipole 
expansion to be viable. 

We choose coordinates such that the motion is restricted to the (x^, x^)-plane. In a 
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FIG. 1: Dipole power as a function of the model parameter a for the data (j45p . The value of 



the quadrupole power as predicted by GR, = 2.248 • 10"^^ W, is indicated by the dashed line. 



corotating frame, the mass density can be written as 

p{x) = 6{x'^)6{x^) [mi6 (a;^ - d/2) + (a;^ + d/2)] 



The dipole moments are therefore given by 
d{mi — 1712) cos{(jjt) 



M\t) 



M\t) = (t-—) , M=^(t)=0, 
V 2u/ 



leading to 



(41) 



(42) 



(43) 



1 



Taking the average of the dipole contribution over one period T := 27i/u we obtain 



9c 



(44) 



dt 



Ovrc jQ 

In Figure 1 we plot Pd against the model parameter a for the data 



mi = 1.337 Mq, 
CO = 0.7112- 10-S-\ 



ma = 1.250 Mq, 
r = 10^°m, 



(45) 
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which is inspired by the PSR J0737-3039 system, but moving in a circular orbit with angular 



velocity 27r/T, where T = 8.834 ■ 10^ s is the measured orbital period 18|. The parameter 



a lies within the interval X = [2 ■ 10^'' m^, 3.25 ■ 10^^ m^], such that the dipole term is com- 
parable to the value of the quadrupole power as predicted by GR, = 2.248 • 10^^ W. 
From f l40p we infer that d ^ 9.5 ■ 10^ m for a G X; the Newtonian value is d = 8.798 • 10^ m. 
We choose the distance r from the source such that d/r ^ 1/10. For this data, the dipole 
power equals 1% of the GR value at a ^ 1.7 ■ 10^^ m^. 



VII. CONCLUSION 

We have derived the f{R) correction terms to the GR quadrupole formula for the emission 
of gravitational radiation to leading order. An important result is that, in contrast to GR, 
quadratic f{R) theory predicts the radiation of monopoles and dipoles. This is the case for 
nearly every alternative metric gravity theory known today [l|. 

When considering a hypothetical binary similar to the double pulsar PSR J0737-3039, 
we found for the constant a an upper limit of about 1.7 ■ 10^'' m^. It would be interesting to 
study the result for real data of pulsars to see whether this limit can be improved. Notice 
that the limit on a we got from the geodetic precession using the Gravity Probe B data 
and the precession of the pulsar PSR J0737-3039B is somewhat more stringent pJi]. When 



completing our paper we came aware of the recent publication 9| as well as the preprint j20| , 



which contain equivalent results to some extent. The multipole expansion we focus on in 



our paper is, however, complementary to 9|], and the treatment in 20] is based on different 
assumptions concerning the conservation laws. 
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